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I. INTRODUCTION 



This paper reports results from the first in a series of planned/proposed experiments to 
demonstrate the formation of imploding spherical plasma liners via an array of merging 
high-Mach-number (M) plasma jets. Results are obtained on the Plasma Liner Experi- 
ment (PLX), depicted in Fig. [U at Los Alamos National Laboratory. Imploding spherical 
plasma liners have been proposed^- as a standoff compression driver for magneto-inertial 
fusion (MIF) 4 ~ 6 and, in the case of targetless implosions, for generating cm-, [is-, and Mbar- 
scale plasmas for high energy density (HED) physics^ research. Several recent theoretical 
and computational studies have investigated the physics of imploding spherical plasma liner 
formation^— and also the fusion energy gain of plasma liner driven MIF.-^^^^ In this 
work, we provide a detailed experimental characterization of high-M argon plasma jet prop- 
agation. In a forthcoming paper we will present the experimental characterization of two 
such jets merging at an oblique angle. The next step, a thirty-jet experiment to form spher- 
ically imploding plasma liners that reach 0.1-1 Mbar of peak pressure at stagnation, has 
been designee?^ but not yet fielded. The primary objective of the single-jet propagation 
and two-jet oblique merging studies is to obtain critical experimental data in order to (1) un- 
cover important unforeseen issues and (2) provide inputs to and constraints on numerical 
modeling efforts aimed at developing predictive capability for the performance of imploding 
spherical plasma liners. 

This study focuses on addressing issues for the potential use of railgun-driven high-M 
plasma jets for forming imploding spherical plasma liners and the ability to reach HED- 
relevant stagnation pressures (> 1 Mbar). Thus, we are interested not only in the jet 
parameters as they exit the railgun, but also in the evolution of the plasma jet as it prop- 
agates a distance of ~ 0.5 m, as this evolution will affect subsequent jet merging and 
ultimately the plasma liner formation and implosion processes. As such, our work con- 
stitutes a unique contribution to the large body of railgun research, 17 which has primarily 
focused on the dynamics and performance of the "armature"— ~— within the railgun bore and 
the ability of railguns to launch solid projectiles or develop thrust for military^ and space 
applications.— ~— These jets, if injected into a tokamak plasma, may also find applications 
in core re-fueling or edge-localized-mode (ELM) pacing.— 

The remainder of the paper is organized as follows: Sec. [Ill gives a brief overview of the 
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physical steps of imploding spherical plasma liner formation using merging plasma jets, and 
summarizes key recognized issues (focusing on jet propagation prior to jet merging); Sec. HTT1 
describes the experimental setup and diagnostics; Sec. IIVI presents the experimental results 
on plasma jet parameters, structure, and evolution; and Sec. [V] provides conclusions and a 
summary. 

II. USE OF PLASMA JETS FOR FORMING IMPLODING SPHERICAL 
PLASMA LINERS 

A. Steps of plasma liner formation 

Here, we provide a brief summary of the steps in imploding plasma liner formation using 
an array of merging plasma jets. Much more detailed accounts, including both theoretical 
and computational results, have been presented elsewhere .-^^ 

First, multiple plasma jets are launched radially inward from the periphery of a large 
spherical vacuum chamber. The jets propagate separately until they coalesce at the merging 
radius i?, m , which depends on the jet M, the number of jets N, the initial jet radius r^o, and 
the chamber radius R w . For non- varying M and the assumption of a jet radial expansion 
speed of 2C s /(7 — 1),— where C s is the ion sound speed and 7 is the polytropic index, R m 
has been derived as^ 

rjo (Mg + 1) + fl m 
"~- 1 + ^(^ + 1)' W 

Note that as M — > 00 (i.e., no radial jet expansion), R m — > N 1 / 2 rj /2. For PLX-relevant 
values (M = 14, iV = 30, r\, = 2.5 cm, R w = 111 cm, 7 = 1.4), R m w 50 cm, meaning that 
the jets will propagate about 60 cm before merging with adjacent jets. The reduced value of 
7 = 1.4 (below the ideal gas value of 5/3) is used throughout the paper to approximate the 
internal degrees of freedom of an argon plasma, i.e., due to ionizations and excitations.— 1 ^ 
If M increases during jet propagation due to radiative cooling, then Eq. (GQ) overestimates 
R m and underestimates the jet propagation distance before merging. Nevertheless, for PLX, 
we are interested in jet propagation distances of order 0.5 m. Characterizing the evolution 
of jet parameters over that distance is the focus of this paper. 

Adjacent jets merge at oblique angles at R m to form an imploding spherical plasma liner. 
For PLX, the jet merging angles are determined by the vacuum chamber port positions, with 
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nearest neighbor jets meeting at 9 — 24°. Even at this angle, two adjacent jets, each with 
M = 14, meet with a relative Mach number of 2Msin(#/2) ~ 6. Thus, shock formation and 
associated heating may be expected to contribute significant non-uniformities to the plasma 
liner formation process. The strength of the shock, and indeed whether shocks even form, 
are open research questions due to the parameter regime of our jets,->22 i.e., the plasma 
within each jet is highly collisional (argon ion mean free path Aj ~ 4 x 10 -4 cm compared 
to the jet diameter of 5-20 cm, where rii = 2 x 10 16 cm -3 , Tj = 1.4 eV, and mean charge 
state Z e ff = 1 have been used to estimate A«) but the interaction between two jets is semi- 
collisional or even collisionless (due to the high relative velocity between the jets) on the 
scale of the jet diameter. Shock heating, should it occur, may degrade the implosion Mach 
number, which would lead to the undesirable result of lower liner stagnation pressure.— 'i 3 - 
The non-uniformities introduced by discrete jet merging may lead to asymmetries in the 
plasma liner implosion that significantly degrade the peak achievable stagnation pressures. 
These and other issues relating to plasma liner formation via discrete merging plasma jets 
have been studied theoretically 8 and computationally 11 elsewhere and are beyond the scope 
of this paper. Experimental results from PLX on two-jet oblique merging will be presented 
in a forthcoming paper. 

Finally, after all the jets merge into an imploding spherical plasma liner, the liner con- 
verges with mass density p ~ poR~ 2 , where po is the mass density at R m and R is the radial 
position of the imploding liner.-* 1 ^^ This results in the convergent amplification of liner 
ram pressure until the liner reaches the origin and stagnates (for targetless liner implosions), 
converting the liner kinetic energy into stagnation thermal energy and excitation/ionization 
energy of the liner plasma. 13 An outward-propagating shock is launched, and when this 
shock meets the trailing edge of the incoming liner, the entire system disassembles after a 
stagnation time T stag ~ Ljo/VjQ,—^- where Lj and Vjo are the initial jet length and velocity, 
respectively. 

B. Issues relating to jet propagation 

There are several issues relating to jet propagation that are important for determining the 
performance of the subsequent liner formation, implosion, and stagnation. Assessing these 
issues experimentally is the primary motivation for this single-jet study. In all cases, the 
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experimental results are intended to improve the predictive capability of imploding spherical 
plasma liner modeling. 

The first issue is the achievement of the requisite jet parameters for reaching the desired 
liner stagnation pressure. The design goal for PLX is to reach 0.1-1 Mbar with a total 
liner kinetic energy of about 375 kJ. A 3D ideal hydrodynamic simulation study™ exploring 
a wide parameter space in plasma jet initial conditions contributed to the PLX reference 
design calling for thirty argon plasma jets with initial density ~ 10 17 cm -3 , velocity w 
50 km/s, and mass ~ 8 mg. The simultaneous achievement of these parameters has recently 
been demonstrated at HyperV Technologies, to be reported elsewhere. For this paper, we 
operated at reduced values to extend the lifetime of the railgun and maximize the number 
of shots. The second issue is the evolution of jet parameters, especially density and velocity, 
during propagation. Developing an accurate predictive capability for modeling plasma liner 
formation via merging jets requires accurate knowledge of these jet parameters, both at the 
exit of the railgun as well as at R m . The jet velocity is expected to remain nearly constant, 
but density decay arises from jet radial and axial expansion. The third issue is the nature 
of the jet radial and axial profiles because they are important for accurate modeling of the 
jet merging and liner formation/implosion processes. 

Experimental values (presented in Sec. IIVP allow the jet expansion to be estimated. In 
a purely hydrodynamic treatment in which the jet has a non-varying M, both the radial 
and axial expansion speeds of the jet can be estimated to be between C s (jet bulk) and 
2C s /(7 — 1) (jet edges).— For argon at T e = 1.4 eV and 7 = 1.4, these values are 2.2 km/s 
and 11 km/s, respectively. If the jet travels at 30 km/s over a distance of 50 cm, this gives 
a transit time of 16.7 /is, during which the jet radius will increase by 3.7-16.7 cm and the 
jet length by 7.4-33.4 cm. For a jet with rjo = 2.5 cm and Ljq = 20 cm, the jet volume 
will increase by a factor of 8.4-157.5, and the density will drop by the same factor. Due 
to the nearly factor of twenty uncertainty in the theoretically predicted density decay, and 
the fact that jet cooling due to expansion and radiative losses (and thus a varying rate of 
expansion) are not accounted for in the above treatment, it is imperative to determine the 
density decay by direct experimental measurement. 
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III. EXPERIMENTAL SETUP 



A. Plasma Liner Experiment (PLX) 

Experiments on PLX are conducted in a 9 ft. (2.74 m) diameter stainless steel spherical 
vacuum chamber, as depicted in Fig. [TJa), situated in a 3000 ft. 2 (279 m 2 ) high-bay space 
with a ten-ton overhead crane. The vacuum chamber has 60 smaller ports (11 in. outer and 
7| in. inner diameters) and 10 larger ports (29.5 in. outer and 23.5 in. inner diameters); all 
flanges are aluminum. We presently have two operational plasma railguns (see Sec. IIII Bl for 
railgun specifications) installed on the vacuum chamber to study high-M single-jet propaga- 
tion, two-jet oblique merging, and two-jet head-on merging. This paper reports only single 
jet results. The vacuum base pressure is typically in the low-10~ 6 Torr range, achieved 
with a turbo- molecular pump (3200 1/s Leybold Mag W 3200C, 12 1 in. pumping diame- 
ter) backed by an oil-free mechanical pump (Edwards IQDP 80). During each experimental 
shot, the turbo-pump is isolated by closing a gate valve, and thus the vacuum pressure 
is in the 5-30 /iTorr range when a plasma jet is fired into the chamber. The pressure in 
the chamber after a shot is in the 0.3-to-few mTorr range; the gate valve is opened and 
the chamber evacuated prior to the next shot. All chamber pressures are recorded using 
an MKS 972B DualMag transducer. A LabVIEW-based (http://www.ni.com/labview) 
40 MHz field programmable gate array (FPGA) system controls safety interlocks, shot se- 
quence including bank charging and dumping, and all trigger signals. Sensitive control 
and data acquisition electronics reside within an electromagnetically shielded cage or other 
shielded racks. All time-series data shown in this paper were digitized at 40 MHz sampling 
rate, 12-bit dynamic range, and 1 mV bit resolution (by Joerger model TR digitizers). All 
times reported in the paper are relative to the trigger time of the gun rails. Experimental 
data and shot data are immediately stored into an MDSplus (http://www.mdsplus.org) 
database after every shot. 

B. Plasma railgun 

To achieve the design objectives of PLX (i.e., 0.1-1 Mbar of peak liner stagnation pres- 
sure using thirty plasma jets with total implosion kinetic energy of ~ 375 kJ) within bud- 
getary constraints, two-stage parallel-plate railguns (see Fig. [2]), designed and fabricated 
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by HyperV Technologies Corp.,— were selected as the plasma gun source for forming and 
launching jets with the requisite parameters (jet density ~ 10 17 cm -3 , velocity ~ 50 km/s, 
and mass ~ 8 mg). The development, optimization, and performance scaling results of 
these PLX railguns, as well as the simultaneous experimental achievement of the afore- 
mentioned jet parameters, will be reported elsewhere. Note that larger coaxial guns with 
shaped electro des 3 ^ - — are also being developed due to their suitability for very high current 
(> 1 MA) and high jet velocity (> 100 km/s) operation, and for having attributes that 
minimize impurities, as potentially required for the MIF application. 

The railgun bore cross-sectional area is 2.54 x 2.54 cm 2 . The gun has a fast gas-puff valve 
(GV), a pre-ionizer (PI), HD-17 (tungsten alloy) rails housed in a Noryl (blend of polypheny- 
lene oxide and polystyrene) clamshell body, zirconium toughened alumina (ZTA) insulators 
between the rails, and a cylindrical acrylic nozzle (5 cm diameter and 19 cm length). The 
results reported here are from shots with an underdamped, ringing, and slowly decaying 
gun current (Fig. [3]) that produces multiple jet structures, with the leading structure having 
a length of order 20 cm. We use ultra-high-purity argon (typically at 18-20 psig) for the 
gas injection, injecting an average total mass of about 35 mg in each shot for the data set 
presented in this paper (but only a small fraction of the total injected mass is contained in 
the leading jet structure). 

The railgun firing sequence is GV followed by the PI and finally the rails. The time 
delays are all adjustable and are chosen for optimizing certain aspects of jet performance 
such as density, velocity, or mass. For the shots reported here, the GV is fired 300 /is before 
the gun rails so that neutral argon fills the PI volume and the very rear of the railgun bore. 
The PI is fired 30 fis before the gun rails to break down the neutral gas, giving the PI 
plasma just enough time to fill the very rear of the railgun bore. Then the rails are fired to 
accelerate the PI plasma down the bore. Control over the plasma jet density and mass is 
mainly through the feed line pressure, GV bank voltage, and timing of the neutral gas puff; 
control over the jet velocity is mainly through the railgun current and charge voltage. For 
single-gun operation, the gun, PI, and GV are driven by 36 /xF, 6 fiF, and 24 /iF capacitor 
banks, respectively, charged typically to -24 kV, 20 kV, and 8 kV, respectively. The 36 and 
6 /iF banks use 60-kV, 6-/iF Maxwell model 32184 capacitors, and the 24-/iF bank uses a 
single 50-kV Maxwell model 32567 capacitor. All banks are switched by spark-gap switches 
triggered via optical fibers. Details about the gun design, operation, performance scaling, 
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and best achieved parameters will be reported elsewhere. 

We have evidence that the plasma jet is not 100% argon. We observe an approximately 
25% higher pressure rise in the chamber (for the data set presented in this paper) when the 
GV and railgun are both fired, compared to when only the GV is fired. The pressure discrep- 
ancy is most likely explained by the railgun current ablating material from the HD-17 rails 
and possibly also the ZTA insulators. We also observe hydrogen, oxygen, and aluminum im- 
purity spectral lines, and potentially others that have not yet been identified (with tungsten, 
nickel, iron, and copper from the rails as likely candidates). In the remainder of the paper, 
we have assumed that the plasma jets are 100% argon for the purpose of interpreting the 
experimental data because the error introduced is generally small compared to diagnostic 
measurement uncertainties. Impurity control in plasma jets is clearly an important issue for 
the MIF standoff driver application and requires further study. 



C. Diagnostics 

Plasma jet diagnostics include an eight-chord interferometer, a visible and near-infrared 
(IR) survey spectrometer, an array of three photodiode detectors, and an intensified CCD 
(charge-coupled device) imaging camera. Details of each plasma jet diagnostic system are 
given below. Figure H] shows the diagnostic views in relation to the plasma jet propagation 
path. A discussion of the entire planned PLX diagnostic suite is described in more detail 
elsewhere.™ We emphasize that all plasma jet diagnostic measurements are averaged quanti- 
ties over their viewing-chords, and therefore we report mostly sight-line-averaged quantities. 

Diagnostics for the railgun and pulsed-power systems include Rogowski coils (for moni- 
toring railgun, GV, and PI discharge currents), Pearson current monitors (model 2877) on 
parallel resistors (for monitoring instantaneous capacitor bank voltages), and five magnetic 
probe coils along the railgun bore (for monitoring electrical current propagation down the 
bore). Figure |3] shows representative gun current I gun , gun voltage V gun , and the gun bore 
current hore from the rearmost gun bore magnetic probe. 
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1. Eight-chord interferometer 



An eight-chord fiber-coupled interferometer was designed and constructed for the PLX 
project.— The system uses a 561 nm diode-pumped, solid-state, 320 mW laser (Oxxius 561- 
300-COL-PP-LAS-01079) with a long coherence length (> 10 m). Along with the use of 
single-mode fibers (Thorlabs 460HP) to transport the laser beams to and from the vacuum 
chamber, this allows for the use of one reference chord for all eight probe chords, as long as 
length mismatches between reference and probe chords are much smaller than the coherence 
length. The fiber-coupling allows for relatively simple chord re-arrangements at launch and 
reception optical bread boards mounted on the vacuum chamber. The beam-splitting and 
combination optics on the main optical table do not need to be altered to re-arrange chords. 
Large borosilicate windows are used on both the launch and reception chamber ports. For 
this paper, the eight laser probe beams (3 mm diameter at the plasma jet) were arranged 
transversely to the direction of jet propagation at different distances from the railgun nozzle: 
Z = 35.0-79.5 cm at equal intervals of approximately 6.35 cm (see Fig. |4|). 

The Bragg cell for the interferometer uses a radio frequency generator (IntraAction ME- 
1002) that produces a 110 MHz signal. The output signal from the final photo-receivers 
are passed through bandpass filters (Lark Engineering MC110-55-6AA) at 110 ± 55 MHz 
to filter out higher-frequency components of the heterodyne mixing and lower-frequency 
electrical noise. The filtered signal is decomposed into two signals, I and Q, proportional to 
the sine and cosine of the signal, respectively. The I and Q signals pass through a 40 MHz 
low-pass filter and then are digitally stored (40 MHz, 12-bit resolution, 50 Q input). A 
detailed description of the system design, components, setup, and electronics are reported 
elsewhere.— 

For a partially ionized argon plasma, the interferometer phase shift of our system has 
been derived as 36 



where / = ni/n tot is the ionization fraction, n tot = rii + n n (where rii and n n are the ion and 
neutral densities, respectively), and the integral is over the chord path length. Note that 
A(f) < when / < 0.07235. Figure |5] shows selected contours of constant A0, calculated 
using Eq. (j2J), as a function of / and J n tot dl. 





9 



2. Survey spectrometer 



Spectrally resolved plasma self-emission is recorded with a survey spectrometer system 
consisting of a 5-mm diameter collimating lens (BK-7), a 19-element circular-to-linear silica- 
core fiber bundle (Fiberguide Industries "Superguide G," 10 m long), a 0.275 m spectrom- 
eter (Acton Research Corp. SpectraPro 275) with three selectable gratings (150, 300, and 
600 lines/mm), and a gated 1024-pixel multi-channel-plate array (EG&G Pare 1420). All 
measurements reported here were taken with the 600 lines/mm grating. We have taken 
measurements from about 300-900 nm, but all the data reported in this paper are between 
430-520 nm. Two corrections are applied to the raw spectrum for each shot: (1) attenu- 
ation of the fiber bundle (approximately negative linear slope around -30 dB/km between 
430-520 nm) and (2) pixel-dependent factor associated with the position of the grating in 
the spectrometer. The latter correction amplifies the spectra at low and high pixel-values 
over the raw spectrum by ~ 1.7-5 depending on the exact pixel value. We placed copper 
mesh in front of the collimating lens, as necessary, to keep the peak counts under about 
16,000 to avoid saturating the detector. For the chord position which is Z f=a 41 from the 
railgun nozzle (see Fig. H]), the counts were of order 100, and thus no mesh was used. Before 
each shot day, a spectral lamp was used to record line emission at known wavelengths to 
provide a pixel-to-wavelength calibration. The spectral resolution for the data presented in 
this paper is 0.152 nm/pixel. The diameter of the viewing chord at the position of the jet, as 
imposed by the collimating lens, is approximately 7 cm, which constitutes a nominal spatial 
resolution for the spectroscopy data. The time resolution is 0.45 as as determined by the 
exposure (gate) time of the detector. One spectrum at one time is taken for each shot. 

3. Photodiode array 

A three-channel photodiode array (PD1, PD2, and PD3) is used to collect broadband 
plasma emission for determining jet propagation speed. PD1, PD2, and PD3 collect light 
mostly transverse to the direction of jet propagation at Z = 2.7, 27.7, and 52.7 cm, re- 
spectively (see Fig. Hj). For each channel, light is collected through an adjustable aperture 
positioned in front of a collimating lens (Thorlabs F230SMA-B, 4.43 mm focal length), 
which is connected to a silica fiber that brings the light to the photodiode detectors inside a 
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shielded enclosure. The silicon photodiodes (Thorlabs PDA36A) are amplified with variable 
gain, and have a wavelength range of approximately 300-850 nm. The peak responsivity is 
0.65 A/W at 970 nm. The quoted frequency response decreases with increasing gain. For 
the data reported in this paper, the channels at Z — 2.7, 27.7, and 52.7 cm had gain set- 
tings of 20, 50, and 50 dB, respectively, corresponding to quoted bandwidths of 2.1, 0.1, and 
0.1 MHz. However, we note that the observed rise times for the second and third channels 
are much faster than the quoted 0.1 MHz bandwidth would dictate (see Sec. II V A lj) . All 
three channels had aperture openings of < 1 cm, which constitutes a nominal spatial reso- 
lution for the photodiode data. The continuous (in time) photodiode signals are digitized 
at 40 MHz by the Joerger TR. 



4- Fast-framing CCD camera 

An intensified CCD camera (DiCam Pro ICCD), with spectral sensitivity from the UV to 
near-IR, is used to capture visible images of the plasma jet. The camera records 1280 x 1024 
pixel images with 12-bit dynamic range. The camera can record up to two images per 
shot. However, for inter-frame times below several tens of /is, the second image often has 
"ghosting" from the first frame. Thus, in this paper, we recorded only one frame per 
shot. The exposure (gate) is 20 ns. The camera is housed inside a metal shielding box 
and mounted next to a large rectangular borosilicate window on the vacuum chamber. A 
zoom lens (Sigma 70-300 mm 1:4-5.6) was used for recording the images shown in this 
paper. The camera is triggered remotely via an optical fiber. To infer quantitative spatial 
information from the CCD images, we recorded images with a meter stick held to the end 
of the railgun nozzle. Using these calibration images, we are able to determine pixel-to- 
centimeter conversion formulae (see Sec. IIVB2j) . The CCD images presented in this paper 



are shown on a logarithmic scale in false color. 



IV. EXPERIMENTAL RESULTS 

In this section, we present experimental results on the jet parameters (including velocity, 
density, temperature, and ionization fraction), structure, and evolution for the experimental 
data set spanning PLX shots 737-819. 
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A. Jet parameters 



1 . Velocity 

The plasma jet velocity Vj e t is determined from photodiode array data and corroborated 
with interferometer data. The velocity is calculated by dividing the distance between viewing 
chords by the difference in arrival times of the peak signal. We use the arrival time of the 
peak signal rather than the leading edge to obtain a more robust estimate of the bulk jet 
velocity. 

Figure [7] shows the three photodiode signals for a representative shot (744) from the data 
set analyzed in this paper. The signals are collected along viewing chords intersecting the 
jet propagation axis at distances Z = 2.7, 27.7, and 52.7 cm, respectively, from the end of 
the gun nozzle (see Fig. H]). For this shot, the peak arrival times are 22.9, 31.6, and 39.8 us, 
corresponding to average velocities of V12 = 28.7 and V23 = 30.5 km/s between the first and 
second pairs of photodiode viewing chords, respectively. For the entire data set considered 
in this paper, V12 = 28.9 ± 3.9 km/s and V23 = 29.4 ± 4.5 km/s, where the uncertainty is 
the standard deviation over the data set. The velocities did not vary significantly across 
this data set due to the relatively narrow range of peak gun current (255-285 kA) and total 
mass injected into the chamber (pressure rise 1.35-1.75 mTorr, corresponding to 31-40 mg 
of argon). A similar analysis using the interferometry data (see Fig. [8](a) for an example) 
from the chords at Z = 35.0 and 47.7 cm yields an average jet velocity of 34.8 ± 5.6 km/s 
(for the entire data set), which compares well with V23 = 29.4 ± 4.5 km/s given above. The 
jet velocity appears to be nearly constant as it propagates over a nearly 50 cm distance. 
At 30 km/s, the jet M = 14 (assuming that the jet is 100% argon, 7 = 1.4, Z e ff = 1, and 
T e = 1.4 eV). 

2. Density, temperature, and ionization fraction 

Plasma jet density n, electron temperature T e , and ionization fraction / are determined 
via a combination of experimental measurements and interpretation with the aid of theo- 
retical analysis and atomic modeling. From spectroscopy, we determine electron density n e 
via Stark broadening of the impurity hydrogen Hg line (486.1 nm), and we estimate T e by 
comparing measured and calculated non-local-thermodynamic-equilibrium (non-LTE) argon 
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spectra. Steady-state, collisional-radiation calculations with single-temperature Maxwellian 
electron distributions in the optically thin limit were performed using the PrismSPECT 
code with DCA (direct configuration accounting). 37 The following processes were included: 
electron-impact ionization, recombination, excitation, de-excitation, radiative recombina- 
tion, spontaneous decay, dielectronic recombination, autoionization, and electron capture. 
The atomic model for argon consisted of 16,000 levels from all ionization stages with ap- 
proximately 1500 over the lowest three ionization stages. The non-LTE calculations also 
provide the ionization fraction / as a function of n e and T e . The / value is then used, in 
conjunction with Eq. (J2J) and the interferometer data, to provide an independent determina- 
tion of n to t = n-i/f- If we assume that n e = rij (i.e., all ions are singly ionized, a reasonable 
assumption at our temperatures and densities), then n e = rij = fn tot . 

Stark broadening of the line is analyzed for the nozzle and chord views (see Fig.HJ). The 
latter coincides with the position of the Z — 41.4 cm interferometer chord. Figure [9] shows 
examples, from two shots (nozzle and chord views, respectively), of how n e is determined 
by fitting the convolution of the measured instrumental broadening (point spread function) 
and a Lorentzian profile to the experimentally measured Hg spectral feature. Curve-fitting 
is performed using the IDL (Interactive Data Language, http://www.exelisvis.com/IDL) 
routine curvef it. The n P is determined via 



n e (cm" 3 ) = 2.53 x 10 14 



FWHM (nm) 

Oil/2 



3/2 



1.50 x 10 13 



FWHM (pixels) 

«l/2 



3/2 



(3) 



where FWHM is the full- width half- maximum of the Lorentzian fit, and a\/2 (in our 
case = 0.085) is the so-called reduced half-width that scales as line shape and has been 
tabulated 3 ^ for many hydrogen lines for the temperature range 0.5-4 eV and density range 
10 14 -10 18 cm -3 . The value of 0.152 nm/pixel for our spectrometer system has been used in 
Eq. 03]). Table H] summarizes our n e results obtained by Stark broadening analysis, and their 
corresponding viewing chords, times, and relative positions in the jet. The n e m 2 x 10 16 cm -3 
near the gun nozzle falls by approximately one order of magnitude after the jet propagates 
approximately 41 cm. Unfortunately, the number of shots on which we could successfully 
perform the Stark broadening analysis was limited because the Hg line typically is obscured 
by a stronger nearby argon line. In the shots analyzed, the nearby argon line is subtracted 
out prior to performing Stark broadening analysis on the H^ line. 

Next, we use the spectroscopy results to obtain an estimate of T e . Figure [10] shows 
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spectrometer data for the (a) nozzle position at t = 17 fis (corresponding to the rising edge 
of the jet) and (b) chord position at t = 36 /is (corresponding to the jet bulk). The figure 
also shows calculated PrismSPECT non-LTE argon spectra for comparison. The Ar II lines 
that appear in the experimental data only appear for T e > 1.4 eV in the PrismSPECT 
non-LTE calculations, thus providing a lower bound estimate of the peak T e in the jet. 
The calculations show that / = 0.96 (for n e = 1.5 x 10 16 cm -3 and T e = 1.4 eV) and 
/ = 0.94 (for n e = 2 x 10 15 cm -3 and T e = 1.4 eV) for the nozzle and chord positions, 
respectively. An upper bound on T e is also estimated based on the argon spectra centered 
around 794 nm (not shown). An Ar II line visible in the calculations at n e = 2 x 10 16 cm -3 
and T e = 1.7 eV, but not visible in our measurements (shots 569-572 and 655-706), indicates 
an upper bound on the peak T e of 1.7 eV for the nozzle view. This upper bound also 
justifies our assumption that all argon ions are singly ionized. For the rest of the paper, we 
use T e — 1.4 eV, which is a more appropriate jet-averaged value, in estimating T e -dependent 
quantities. Although Tj has not been directly measured yet, it is reasonable for the purposes 
of this paper to assume that Tj « T e because the ion-electron thermal equilibration time 
v~ x — (3-2 x 10~ 9 Z^jj In Ari e //iT 3 / 2 ) -1 , estimated to be 0.2 /is, is very fast compared to the 
jet evolution occurring over many tens of fis (where we have used Z e ff — 1, In A = 4.74, 
n e = 2 x 10 16 cm -3 , fi = mi/m p = 40, and T = 1.4 eV in estimating i> e ). 

Finally, using the chord position / = 0.94 determined above, an estimated interferometer 
chord path length « 10 cm at the Z = 41.4 cm chord (see Sec. HVBl) . and Eq. (j2J), we obtain 
an independent determination of n e = fn to t from the interferometry data. The results 
are shown in Fig. [HJ showing reasonable agreement in n e obtained via spectroscopy and 
interferometry. Table HT1 summarizes all the experimentally measured jet parameters at both 
the nozzle and chord positions. 

3. Remark on the jet magnetic field 

The evolution of the jet magnetic field B after the jet exits the railgun has not been 
measured yet on PLX. The gun bore magnetic probes show that the field strength within the 
gun bore is on the order of several Tesla. The classical diffusion time of B is approximately 
t d ~ Ho3 2 /v±i where fiQ = Arc x 10 7 H/m is the vacuum permeability, 5 ~ 2 cm the 
gradient scale length of B, and rf± the perpendicular Spitzer resistivity. For T e = 1.4 eV, 
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n e = 2 x 10 16 cm" 3 , Z eff = 1, and In A = 4.74, rj± = 2.95 x 10^ 4 fim and td = 1.7 //s. Thus, 
for Vj e t ~ 3 cm/ /is, the magnetic field will decay to e~ 3 = 0.05 of its original value after 5.1 /xs 
or 15.3 cm of propagation. Assuming a value of 3 T in the gun bore, the decayed value after 
5.1 /is will be 0.15 T, and the ratio of jet magnetic energy density B 2 /2fi = 9.0 x 10 3 J/m 3 
to kinetic energy density pV 2 et /2 = 3.01 x 10 5 J/m 3 (assuming an argon density of 10 16 cm -3 
and Vjet = 30 km/s) will be about 0.03. Thus, it is reasonable to ignore (to leading order) 
the effects of the magnetic field on jet evolution over 0.5 m. We defer the direct measurement 
of jet magnetic field evolution to future work. 



B. Jet structure and evolution 

In this sub-section, we present results on jet structure via CCD image, photodiode, and 
interferometer data, including quantitative results on jet length inferred from photodiode 
and interferometer data, and jet diameter inferred from CCD images and interferometer 
data. Note that our jets have a primary leading structure with several trailing structures, 
as seen in the Z = 2.7 cm trace after t = 40 fis in Fig. [7J By triggering a crowbar circuit 
to eliminate the ringing current (results not shown in this paper), we have successfully 
eliminated the trailing jet structures in both the photodiode and interferometer data. In 
this paper, we focus only on the properties of the leading jet structure in non-crowbarred 
shots. 



1. Jet length and axial profile 

We estimate the jet length L using the full-width at 1/e of the maximum of the photodiode 
(see Fig. [7J as an example) and interferometer (see Fig. 13(b) as an example) signals versus 
time to get a Atj et for each photodiode and interferometer signal. It follows that L « 
Vj e t At j et , where Vj e t is determined by the difference in arrival times of the signal peaks, as 
described in Sec. IIV A 11 The photodiodes give L\ = 20.6 ± 3.5 cm, L 2 = 41.8 ± 4.4 cm, 
and L 3 = 48.9 ± 7.5 cm (where the variation is the standard deviation over the data set) 
for the photodiode views at Z = 2.7, 27.7, and 52.7 cm, respectively. The velocities used 
are V 12 = 28.9 km/s, (V l2 + V 23 )/2 = 29.1 km/s, and V 23 = 29 A km/s (see Sec. HVAip for 



calculating Li, L 2 , and L 3 , respectively. The interferometer chord at 41.4 cm, which is in 
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between the photodiode views at Z = 25 and 50 cm, gives L int = 46.7 ± 6.5 cm, consistent 
with L 2 < L int < L 3 , where the velocity used is 34.8 km/s. 

The rate of jet length expansion L is determined from both the photodiode and interferom- 
eter data by taking AL/At between adjacent measurement positions. From the photodiode 
data over the entire data set, L\2 = 24.4 ± 5.4 km/s and L23 = 9.0 ± 7.3 km/s, which corre- 
spond to the spatial ranges Z = 2.7-27.7 cm and Z = 27.7-52.7 cm, respectively. From the 
interferometer data, Li nt = 20.2±14.0 km/s for the spatial range Z = 35.0-41.4 cm. If we as- 
sume L/2 = 2C s j (7 — 1) from hydrodynamic theory— where C s = 9.79 x 10 5 (7Z e //T e //i) 1 / 2 , 
we get an independent estimate of T e = 1.8 eV (for L = L 12 = 24.4 km/s, /1 = 40, Z e fj = 1, 
and 7 = 1.4) that is in reasonable agreement with the range of T e = 1.4-1.7 eV determined 
in Sec. IIVA2I 

The jet axial profile can be inferred, qualitatively, from the photodiode and interferom- 
eter signals versus time. The time dimension can be converted approximately to a spatial 
dimension by multiplying by Vj e t ~ 30 km/s or 3 cm//is (which was done above to obtain 
L = Vj et At j et ). The qualitative jet axial profile as it comes out of the railgun is best observed 
in the Z = 2.1 cm photodiode trace (Fig. [7]), which shows a sharp rising edge lasting about 
2 \is (6 cm). This is followed by some structure over about 5 [is (15 cm) and then by a falling 
tail that decays to about 10% of the peak intensity over about 10 /xs (30 cm). Figure [TT] 
gives a good representation of the jet axial profile farther downstream [Z = 41 cm), based 
on interferometer density data. 

For an average Vj et = V\i = 28.9 km/s (see Sec. IIV A lj) as the jet exits the nozzle, this 
corresponds to an average Atj e t ~ -^1/^12 = 7.1 /is, which corresponds to the half-period of 
I gun (see Fig. [3]). This suggests that the jet length is inversely proportional to the frequency 
of the pulsed-power railgun electrical circuit, and that the jet length can be tailored via 
electrical circuit parameters. 

2. Jet diameter and radial profile 

The jet diameter D is estimated by two methods. The first method is examining CCD 
image line-outs perpendicular to the jet propagation direction Z; this method also provides 
information regarding the jet radial profile. The second method uses the decay of the peak 
phase shift A<f) pea k of the eight interferometer chords and the assumption of conservation 
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of total jet mass to deduce the jet diameter. We recognize that the two methods provide 
related but fundamentally different information, i.e., the first method relies on the intensity 
of emission from the jet plasma (which depends on both density and temperature) while the 
second method derives from the interferometer phase shift (which depends predominantly 
on the line-integrated plasma density). Nevertheless, the two methods give us meaningful 
quantitative estimates of the jet diameter and also provide qualitative information on the 
shape of the radial profile. 

Figure [T27 a) shows an example of a CCD image vertical line-out at Z = 41. 4 cm, and the 
estimate of D based on the full- width at 1/e of the maximum. The time of this CCD image 
(t = 36.0 /xs) corresponds approximately to when the peak emission reaches Z « 41 cm 
(based on photodiode data). The line-out is taken after the CCD image is rotated slightly 
(« 2°) such that the jet propagation axis is horizontal. The horizontal-pixel value (x) 
of the CCD image is converted to the Z coordinate (cm) based on the formula Z(x) = 
68.5 — 10 5 x 2 — 0.041x (where < x < 1023). In addition, the cm/pixel conversion AY in 
the vertical direction of the CCD image is obtained using the formula AY(x) = 2 x 10~ 5 x + 
0.041 cm/pixel. The dependence of AY on x is due to the camera perspective, i.e., nearer 
objects appear larger in the image. These formulae were determined from a CCD image of 
a meter stick held to the end of the railgun nozzle. Figure fl27 b) shows D versus time at 
Z = 41.4 cm determined using the full-width at 1/e method from a series of CCD images 
(shots 780-819). Note that D determined via this method is fairly constant, with D ~ 10 cm 
from t = 26-38 /xs, but then D increases after 38 /is, which corresponds to the time when 
the peak of the leading jet structure is passing through Z = 41 cm (which can be deduced 
from Fig. [7]). The expansion in D at later times implies a lower M, suggesting that the 
trailing part of the jet is either hotter, slower, or both. 

We also consider CCD image vertical line-outs at different Z positions from a single shot 
at a single time (t = 30.0 /xs), as shown in Fig. [T2T c). The diameters as determined by the 
full- width at 1/e of each peak are 8.2, 8.5, 8.5, 9.1, and 10.6 cm for Z = 5, 10, 15, 20, and 
25 cm, respectively. These emission line-outs provide information on the radial profile of the 
jet, and can provide further information on n e and T e profiles if compared with synthetic 
emission profiles generated using spectral modeling codes such as Spect3D. 39 

Next, we evaluate jet radial expansion using interferometer data. Figure Ufa) shows that 
A(p pea k decreases with increasing Z chord position. Figure Wjo) shows A<p peak versus time 
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for each chord, and fits a quadratic function in time to the data points. From Eq. ([2]), it is 
apparent that A(f> pea k ~ n tot D by assuming J n to tdl = n tot D and constant /. By invoking 
conservation of total jet mass, we obtain the relationship 36 



n tot {t)D{t) 2 L(t) ~ A<f) peak {t)D{t)L{t) = constant D(t) oc — , (4) 

A0 peafc (t)L(t) 



giving -D(t) in terms of experimentally measured quantities A0 pea fc(t) and L(t), which are 
both shown in Fig. E](b) for shot 744 as an example. By using the analytic fits to A<j) pea k(t) 
and L(t) given in the legend of Fig. |8](b), it is straightforward to calculate D(t) to within a 
constant that can be determined using D(t = 36 /j,s) ~ 10 cm from Fig. [121(b) . The result 
is shown in Fig. [T5J Note that the diameter obtained at Z m 41 cm obtained from CCD 
image line-outs is likely an underestimate because the jet emission falls off more sharply than 
density if there is a peaked temperature profile. We regard the uncertainty of our reported 
diameter results based on CCD line-outs to be around a factor of two. The nominal radial 
expansion rate is ~ 7 km/s (between 36 and 43 /is), which is within the range L/2 = 
10.1 ± 7.0 km/s obtained between the interferometer chords at 35.0 and 41.4 cm. 

Finally, we show examples of the plasma jet radial density profiles deduced using Abel 
inversion^ of two chord-arrays of interferometer data from a single shot (1106). This shot is 
not part of the data set analyzed in the rest of this paper because the Abel inversion analysis 
requires that the interferometer chords be arranged with different impact factors relative to 
the jet propagation axis Z, as shown in Fig. IT4T a). For the main data set in this paper, the 
interferometer chords were arranged to intersect the Z axis at different values of Z. Our 
Abel inversion analysis assumes cylindrical symmetry for the jet and four radially concentric 
zones of uniform density in each zone. The radii of the center of the zones correspond to the 
four interferometer chord positions (impact factors), respectively, for each Z location, i.e., 
5, 10, 15, and 20 cm for the Z ps 60 cm chord array. The results of the Abel inversion are 
shown in Figs. [14Tb) and IT4Tc). which show (/ — fo)n to t and radial line-outs of (/ — fo)n to t 
and n e , respectively. From Fig. [T4T c). it can be seen that the jet density profile is about a 
factor of two larger than the density profile from CCD image line-outs shown in Figs. [127 a) 
and (He). 
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V. CONCLUSIONS AND SUMMARY 



The issues we address in this work (Sec. Ill Bp are determining the jet parameters, eval- 



uating the evolution of the jet parameters as the jet propagates over about 0.5 m, and 
characterizing the jet axial and radial profiles to the extent possible. These issues are im- 
portant for accurate assessments of the formation of imploding spherical plasma liners using 
merging plasma jets. The experimental data, while not in all cases definitive due to diag- 
nostic limitations, offer crucial information about jet initial conditions and constraints on 
jet evolution that can enhance the accuracy of numerical modeling predictions.— MdLM We 
reiterate that our plasma jet diagnostic measurements are all averaged over their viewing- 
chords, and therefore our reported results are mostly chord-averaged quantities. 

We have presented results on jet parameters, summarized in Tables HI and Hit showing that 
we are within a factor of 2-5 of what is needed to field the thirty-jet imploding plasma liner 
formation experiments specified in the PLX design.->iA Experiments at HyperV Technologies, 
using the same railgun design as the railgun used in this work but operating at > 500 kA, 
have demonstrated the simultaneous achievement of the full PLX design parameters: density 

10 17 cm" 3 , velocity ~ 50 km/s, and jet mass ~ 8 mg. Those results will be reported 
elsewhere. We operated at reduced railgun current (< 300 kA) in order to extend the 
railgun lifetime and maximize the amount of data collected. Although not emphasized 
in this paper, we have also learned how to operate the railgun reliably and taken note of 
where technological improvements should be made (mostly in the details of the pulsed-power 
system) for a thirty-gun experiment. 

The need to understand the evolution of the jet as it travels a distance of about 0.5 m 
was a primary motivation for this work. To this end, our diagnostics were focused on two 
locations: the nozzle position right at the nozzle exit and the chord position about Z = 41 cm 
away from the nozzle. The diagnostic views were transverse or mostly transverse to the jet 
propagation axis. We were able to determine that average n e fell by about one order of 
magnitude (from 2 x 10 16 to 2 x 10 15 cm -3 ) over this distance, while average T e and / remained 
nearly constant (within our measurement resolution). We were also able to determine the 
approximate jet length and diameter (see Table Hill . and their evolution, showing that the jet 
volume increased by about one order of magnitude which is consistent with the drop in n e . 
As discussed in Sec. IHBl the purely hydrodynamic prediction of constant-M jet expansion 
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has a large uncertainty due to a lack of precise knowledge of the exact expansion rate. This 
uncertainty is further exacerbated by the theory not accounting for radiative cooling which 
is expected to be important in our parameter regime. Using our measured jet parameters, 
the hydrodynamic theory predicts that the jet volume could increase by up to a factor of 
approximately 150. Our experimentally observed volume increase and density drop by about 
a factor of ten provides a useful constraint for validating numerical modeling results (that 
also include the effects of atomic physics) on jet propagation. 

Adiabatic expansion of the jet would dictate that nTV 1 = constant. With a density drop 
of ten, a volume increase of ten, and 7 = 1.4, the average temperature should have dropped 
by a factor of about 2.5. The discrepancy between the latter and what was observed could 
potentially be explained by Ohmic heating associated with magnetic energy dissipation. The 
initial plasma jet thermal /3 is < 0.01, and thus the magnetic energy dissipation (discussed 
in Sec. IIV A 31) could easily balance thermal and radiative losses during jet propagation. 

The radial and axial profiles of the jet are important variables for determining the dynam- 
ics of subsequent jet merging. The profiles are also needed for accurate numerical modeling 
of jet merging and plasma liner formation. We have assessed the radial and axial profiles to 
the extent possible via direct experimental measurements, and they are shown in Figs. [TT1 
fl2T a). fl2T c). and IT4T c). The profiles inferred from CCD image line-outs are dependent on 
both n to t and T e , whereas the profiles inferred from interferometer data are dependent mostly 
on ntot- If T e is peaked in both the radial and axial directions, as can be reasonably ex- 
pected, then this means that the CCD line-outs will underestimate the width of the density 
profile. These profile data provide the opportunity to validate numerical modeling results 
by comparing the experimental profiles with synthetic data from post-processed numerical 
simulation results. 

In summary, we have reported experimental results on the parameters, structure, and 
evolution of high-M argon plasma jets launched by a pulsed-power-driven railgun. An array 
of thirty such jets has been proposed as a way to form imploding spherical plasma liners to 
reach 0.1-1 Mbar of stagnation pressure. Imploding plasma liners have potential applications 
as a standoff driver for MIF and for forming repetitive cm-, /is-, and Mbar-scale plasmas for 
HED scientific studies. 
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shot 


view 


time (/us) position in jet 


n e cm ° 


744 


nozzle 


1 7 




2.2 ± 0.1 x 10 16 


763 


nozzle 


17 


rising edge 


1.6 ± 0.3 x 10 16 


738 


nozzle 


30 


late in jet 


8.6 ±0.0 x 10 15 


771 


nozzle 


35 


late in jet 


5.8 ±0.1 x 10 15 


785 


chord 


28 


rising edge 


2.0 ±0.6 x 10 15 


790 


chord 


38 


just past peak 2.0 ± 0.9 x 10 15 



TABLE I. Summary of sight-line-averaged electron density determined via Stark broadening anal- 
ysis of the spectral line (486.1 nm). Nozzle view corresponds to just outside the railgun nozzle, 
and chord view corresponds to Z ~ 41 cm (see Fig. The chord position results are also shown 
graphically in Fig. [TTJ The quoted uncertainties in n e are the standard deviations obtained from 
curve- fitting. 



nozzle (Z ~ 2 cm) chord (Z ~ 41 cm) 



n e (cm" 3 ) 2 x 10 16 2 x 10 15 

T e (eV) 1.4 1.4 

V (km/s) 30 30 

/ 0.96 0.94 

L (cm) 20 45 

D (cm) 5 10-20 



TABLE II. Summary of plasma jet parameters at the nozzle and chord positions, as determined 
by experimental measurements. The n e and T e values represent jet-averaged values. 
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Vacuum Chamber 




FIG. 1. (a) Schematic of the Plasma Liner Experiment (PLX), designed as a thirty plasma jet 
experiment to form spherically imploding plasma liners. This paper reports single jet characteri- 
zation and propagation studies, (b) Illustration of imploding spherical plasma liner formation via 
the merging of thirty jets, and also the stagnated liner at peak compression (later in time). 
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30 cm 
60 cm 



FIG. 2. (a) Three-dimensional view of the railgun used in this work, (b) Side-view schematic of 
railgun showing (from left to right): fast gas valve, capillary (pre- ionizer), bore (HD-17 tungsten 
alloy rails with zirconium-strengthened-alumina insulators), and acrylic nozzle. The distance from 
the back of the rails to the end of the nozzle is 47 cm. 




FIG. 3. Representative railgun I gun , V gun , and gun bore current hore (near the rear of the rails) 
for the shots analyzed in this paper. 
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Spectrometer 
Chord View 




Spectrometer 
Nozzle View 



FIG. 4. Diagnostic setup for the experiments reported in this paper. The spectrometer chord view 
corresponds to the Z = 41.4 cm interferometer chord, and the spectrometer nozzle view is at the 
exit of the railgun nozzle. PD1, PD2, and PD3 are the three photodiode views. 
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FIG. 5. Contours of constant interferometer phase shift A(j) as a function of ionization fraction / 
and line-integrated argon ion plus neutral density J ntotdl, as calculated using Eq. (|2j). 
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FIG. 6. Plasma jet evolution as recorded by the CCD camera over 8 separate shots (800 and 784- 
790). The railgun nozzle is at the very right edge of each image. The images show the logarithm 
of the CCD intensity in false color. 
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FIG. 7. Photodiode array signals versus time. Horizontal lines denote the duration Atj e t over 
which the signals are greater than 1/e of the peak value (used in calculating jet length and rate of 
axial expansion). In this case, the velocities are 28.7 and 30.5 km/s between the first and second 
pairs of photodiodes, respectively. The Atj e t values are 7.9, 14.9, and 17.4 [is, and the jet lengths 
are 22.8, 44.2, and 53.0 cm at Z = 2.7, 27.7, 52.7 cm, respectively. 
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Phase Shift and Ionization Fraction Dependent Particle Density vs Time 10 « 




35r 



25 



15 



20 30 40 50 60 70 80 90 10 

Time [us] 
(a) 

Measured Expansion Parameters vs Time 



I 


— m- Measured A 6 

^peak 

Polynomial Fit A <> = (2.09e+9)'l 2 - (2.05e+5)'t + 5.17 

— Measured Jet Length 

Linear Fit L(t) = (2.1 1e+6)'t- 36.4 





























75 



36 



38 



40 42 44 

Time [us] 

(b) 



46 



48 



45 



35 



FIG. 8. (a) Interferometer phase shifts versus time for all eight interferometer chords (distances 
Z from the jet nozzle are indicated in the legend). Methodologies for calculating jet velocity and 
length are shown, (b) The peak phase shift A<j) pea k (squares) and jet length (diamonds) are shown 
for each chord at the time at which A<j> pea k for each chord occurs. Also shown are analytic fits to 
each set of data. 
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FIG. 9. Determination of electron density n e via Stark broadening of the Hg line for shots (a) 744 
(t = 17 /US, nozzle view), and (b) 785 (t = 28 [is, chord view). Shown are the experimental 
data (diamonds with error bars = ±\/counts), an overlay of the measured instrumental broadening 
profile (dotted line, labeled as "psf" for point spread function), and a Lorentzian profile (dashed 
line) that gives the best fit (minimum x 2 ) of the convolution (solid line) of the psf and the Lorentzian 
to the data. The n e is calculated from the Lorentzian full-width half-maximum (FWHM) using 
Eq. (|3|), and the uncertainty in n e is based on the standard deviation of the FWHM determined 
from curve-fitting. 
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FIG. 10. Intensity versus wavelength from the survey spectrometer with spectrometer views at 
the (a) nozzle (t = 17 /xs) and (b) chord (t = 36 /xs) positions. Shown in both plots are also 
the argon spectra of non-LTE PrismSPECT calculations with the parameters indicated in the 
respective legends. Appearance of the Ar n lines in the data indicated by asterisks implies that 
peak T e > 1.4 eV according to the PrismSPECT calculations. 
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FIG. 11. Interferometer phase shift A(p averaged over shots 775-819 for the Z = 41.4 cm chord 
(left hand y-axis) and electron density n e (right hand y-axis) versus time. The square data points 
are derived from Eq. ([2]) using / = 0.94, interferometer A<p data from the Z = 41.4 cm chord, and 
the corresponding path length estimate based on the jet diameter D obtained from CCD line-outs 
(discussed in Sec. IIVB 2p . Error bars on the square data points represent uncertainty based on 
doubling the path length, which would halve the average n e . The two diamond data points are 
from Stark broadening analysis of the Hg line from spectroscopy (shots 785 and 790). Each discrete 
data point corresponds to a separate shot. 
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FIG. 12. (a) CCD image line-out versus Y (perpendicular to jet propagation direction) and the 
full-width (FW) at 1/e (dashed line) used to define the jet diameter D; (b) D versus time at 
Z = 41.4 cm from CCD line-outs (shots 780-819); (c) CCD line-outs at different Z positions from 
a single shot, with FW at 1/e jet diameters of 8.2, 8.5, 8.5, 9.1, 10.6 cm for the increasing Z values, 
respectively. 
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FIG. 13. Plasma jet diameter versus (a) time and (b) Z, where the squares are from the data at 
each interferometer chord, and the solid line is from Eq. using the fitting functions for A<p pea k 
and L given in Fig. E^b). Error bars represent ±10% uncertainty in the determination of L, which 
is used in Eq. Q. 
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FIG. 14. (a) Location of interferometer chords overlayed on the CCD image for shot 1106, on which 
we performed an Abel inversion analysis to obtain jet radial density profiles (for the yellow and red 
chord arrays separately), (b) Plots of (/ — fo)ntot [cm -3 ] versus radial chord position (relative to 
the jet axis) and time for the Z ~ 50 cm (top) and Z ~ 60 cm (bottom) chord positions; vertical 
lines indicate the times for which line-outs are shown next, (c) (/ — fo)ntot and n e (data points) 
versus radial chord position at different times; / = 0.94 and a maximum path length of 44 cm 
(corresponding to the path length of the innermost chord at Z « 60 cm) were used in calculating 
n e . Error bars represent uncertainty in the chord positions of ±0.5 cm and also account for the 
slight variation in Z of each chord array. 
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